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Abstract 

Recently, a necessary and sufficient condition to de- 
termine the robust stability of a multilinear interval 
control system has been reported in [1,2] as an exten- 
sion of the well-known Box theorem [3] which deals 
with the linear affine case. This paper introduces a 
simple but computationally efficient algorithm, based 
on the above result, to check the robust stability of 
such systems. The method is also extended to find 
the parametric stability margin of such a system. 

1. ROBUST STABILITY 

Let p := \pi P 2 • * * Pi] a vector of real parameters 
lying in the interval uncertainty set 

n := {p I p- <Pi <pf, i= 1,2, ••■,/}• (1) 

Consider the polynomial 

tf(s,p) := 6o(p)+6i(p)a+Mp)* 2 + ---+Mp)» n ( 2 ) 

wherein the coefficients 6i{p) are affine multilinear 
functions of p. We shall refer to this type of polyno- 
mial as a multilinear interval polynomial. It is easy 
to show that any multilinear interval polynomial can 
be rewritten as 

m r 

*W:=E^WII P «W (») 

.'=1 j - 1 

where Pij(s) are interval and Qi(s) are fixed poly- 
nomials. The necessary and sufficient condition for 
robust stability of the family 6(s) under the assump- 
tion that Pij{s) are independent is given below. Let 
be the set of Kharitonov vertex polynomi- 
als (4] associated with the interval polynomial PiA°) 
and Sp^^yXj) the set of Kharitonov segment poly- 
nomials associated with interval polynomial ^( 5 ). 

Theorem 1. [X ,2] The multilinear family 6(s ) is 
Hurwitz stable if and only if the set of manifolds Mj 
for i = 1, 2, • - • , m are Hurwitz stable where 

m r r 

M/ = ^ Qt(*) IT U ^J') 

*=!,»*! ;=i j = 1 

( 4 ) 


for Xj E [0, 1] for all i. 

Now let us define the new parameter vector 

A •= [Ai ^2 ’ * ’ ^>]» (5) 

then each manifold in M * can be written as the poly- 
nomial set 

6{s } X) := 6 0 {X) + *i(A> + • ■ - , A< € [ 0,1] (6) 

where the coefficients <5* (A) are multilinear functions 
of A- Therefore, the problem of checking robust sta- 
bility of a multilinear interval control system is re- 
duced to checking the stability of the set of multilin- 
ear polynomials shown in eq. (6), as A ranges over the 
positive unit hypercube in the first quadrant. 

Let us define the set 

A := | e [0,1]}. (7) 

Let V denote the vertices of the A set , 

V := {A | Xi = 0 or A< = 1, for all i} (8) 
and Ay := {£(s, A) | A 6 V} (9) 

denotes the set of vertex polynomials. 

Note that A is a continuum of polynomials whereas 
Ay is a discrete set of polynomials and 

Ay C A. (10) 

Fixing s = s*, we let A(s*) denote the set of points 
£(a*,A) in the complex plane obtained by letting A 
range over [0,1]: 

A(0 :={*(**, A) I Ae [0,1]}. (11) 

Likewise we have the discrete set of points in the com- 
plex plane 

A V (0 :={*(•*, A) | AeV}. (12) 

We now state the well known Mapping Theorem [6]. 
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Theorem 2. (Mapping Theorem [0]) Under the 
assumption that 6i (A) are affine multilinear functions 
of A 

co A (s*) = co Ay(s*) (13) 

for each s* 6® where co (•) indicates the convex hull 
of a set (•). 

This theorem shows that the image set of the mul- 
tilinearly parametrized interval family, evaluated at 
any point s* is contained in the convex hull of the 
vertices evaluated at the same point. Although the 
convex hull of the vertices overbounds this set we can 
improve the accuracy of the approximation arbitrarily 
by introducing additional vertices. This is illustrated 
in Figures 1 to 4. The image set A(j*) is contained 
in the convex hull of R(s*) = co Ay (s*) as shown in 
Figures 1 and 2. As shown in the subsequent figures 
(Figures 3, 4) this approximation can be improved by 
decomposing A hypercube as a union of smaller boxes 
and thereby introducing additional vertices as shown 
in Figures A.l and A. 2. 



From this theorem, we can easily establish the follow- 
ing corollary. 

Corollary 1. The multilinear polynomial set A(s) 
is Hurwitz stable if the convex hull o/Ay(i) is Hur- 
witz stable. If the dependence of the coefficients 6»(p) 
on p is linear the stability of A(s) is equivalent to 
that of the convex hull of Ay (s). 

If we define the set of convex combinations of the 
vertex polynomials in Ay(j) by 

E(a) := {j*, •$<(*) + (1 - M», )*;(«) I Mo G [0, 1] and 

M*)» $»'(*) e A v) > ( 14 ) 

the stability of co Ay(s) is equivalent to the stability 
of E(a). Since the set E(s) consists of line segments 
joining every pair of vertices in Ay($), its stability 
can be easily verified by the Segment Lemma [7]. The 
segment lemma basically determines if the phase dif- 
ference of a pair of stable vertex polynomials reaches 


180 degrees at some w. If the phase difference does 
not reach 180° for any w, the line segment joining 
the two vertices is Hurwitz stable. This condition is 
called the Phase Condition. 

Using the above concepts we can easily prove the fol- 
lowing. Let A (ju>) denote the complex plane image 
of the set A(s) evaluated at s = ju>. 

Theorem 3. The set of multilinear polynomials 
A (s) is Hurwitz stable if i) for some w E [0,oo), 
0 £ A (jut) and ii) the set of all corresponding line 
segments E(5) u Hurwitz stable. 

This theorem along with the previously given proce- 
dure for approximating the image set can be used to 
develop an efficient computational technique to check 
robust stability. We can also extend this result to 
the computation of the parametric stability margin 
of the multilinear interval control systems. This is 
done next. 

2. PARAMETRIC STABILITY MARGIN 

Consider a Hurwitz polynomial of the form 

«w = Ea«wn^w ( 15) 

i=l 3=1 

where (?*($) are fixed and 

Pij (a) := Pij + PijS + p^a 2 + ■ ■ ■ . (16) 

Let us assume that coefficients jf-- for k = 0, 1, 2, • • • 
are subject to variations. If we assume the variations 
in coefficients are bounded as 

P?, e ~ w b € ’ Pu„ + w i3 e h £ Z° ( 1? ) 

for a fixed value of the weights Wij with being the 
nominal values of the parameters, then the paramet- 
ric stability margin is defined as the maximum value 
e* so that the multilinear polynomial £(s) remains 
Hurwitz stable for all e E [0,6*). 

Let us recall the manifolds shown in eq. (4). From 
Theorem 1, the parametric stability margin of the 
polynomial in eq. (15) is equivalently defined as 

min {e | max{e | Af* E 7t } VA< E [0, 1]}, 

lEm.iEf}, (18) 

Note that each manifold in Mi has the form of 
Q^K^K^s) • • • K' P \' r (s) + • • • 

(a)xjj- 1 ’ (a) • • • ifc; w 

+g, (5)^(5, x 1 )s i ^(,, x 2 ) ■ • • s% r (»,K) 
+Qi + i(*)Kk;:S')Kiz « • • • + • • ■ 


where K^. (a) is the ib th Kharitonov polynomial and 

Sp^Xa, A) is the fc‘ h segment polynomial associated 
with an interval polynomial P l} (a) , respectively. 

Using the sufficient condition developed in Theo- 
rem 3, we can develop the following computational 
procedure. For some c, if all the corresponding line 
segments in E(s) are Hurwitz stable for all e E [0,e), 
then the actual parametric stability margin e* is al- 
ways greater than or equal to e. However, if we in- 
troduce additional vertex points the approximation 
of the image set by the convex hull of the vertices im- 
proves and so the difference e* — e becomes smaller. 
This permits us to compute the parametric stability 
margin with arbitrary accuracy depending upon how 
much we want to refine the approximation and how 
much computation al burden may be taken. From the 
above considerations, a bisection algorithm is devel- 
oped as follows: 

Bisection Method: e 
For each fixed e 

Find all manifolds for 

(Ai A,) 6 {(0, - - - , 0), (1, - • * , 1)} 

IF all E(s) satisfy phase condition 
THEN forward bisection 
ELSE backward bisection 
ENDIF 

Bisection Method stops when the section piece 
is small enough. 

3. NUMERICAL EXAMPLES 

Example 1. Consider the feedback system shown in 
Figure B. 


-fr 


C(s) 


P(s) 


Figure B. Feedback System. 


diE <£] 
do €[«£,<] 
and , n* ] 

["o » n o 1 
[dj i<^] 
[dj ,d^] 
[do~,do + ] 


with the nominal value = 1-9 
with the nominal value d£ = 2.0 
= [n? - w ni e, n° + u)„, e] 

= [ng -iu no e,n° +u>n„e] 

= [d£ - w d 2 c,(% + to d ,e] 

= [d? - w dl €,Ji + w dl e] 

= [dg -w d <,*,<% + »*>«]• 


Then the characteristic polynomial is given by 


II := {d(s)d*(j) + n(a)n,;(*)} 

which shows that the parameters are entering into its 
coefficient linearly. As we stated in Corollary 1, the 
phase condition becomes necessary and sufficient in 
this case. The sets of Kharitonov polynomials asso- 
ciated with n(s,e) and d(s,e) are defined as follows: 

f ^n(*J C ) = ri O n i Kn{ 3 i € ) ~ n O d" n i 3 

\ K n{*> e ) = "o +«r«. K n{ s > e ) = n O + n l a 

f K}(a,e)= d^+d^a + d+a* 

) K d (a,e) = <£+d+a + d+a 2 

} K d (a,c) = d++d-s + d 2 -s 2 

l^(s,e)= ctf+d+a + dja 2 . 

Thus we need to check the phase condition of the 
following set of vertex pairs: E := Ei U E 2 where 

Ei := {(-Ki(*,c)dr(*) + Ki(a,e)n c (a), 

K^a,e)d r .(a) + Ki(a,€) nn (s))} 

e 2 ~ {(■*:5( a i c K(*) + - K i(*» e ) n «( , )» 

/f5(s,e)d K (s) + «■;(«, e)n fi (a))} 
for i = 1, 2, 3, 4 and 

(t, k) e {(1, 2), (1, 3), (2, 3), (3, 4)} . 

Using a bisection method on e, we have the paramet- 
ric stability margin: 

e* = max{e | E E W, V e E [0, e*]} = 0.146. 


Example 2. Consider the following interconnected 
feedback system shown in Figure C. Let 


n(s) TliS + Tip 

d(s) * d 2 s 2 + dis -f do 

n r (s) s 2 + 2s + 1 

d r (s) ’ s 4 + 2s 3 + 2s 2 + s 

with its plant coefficient variations being bounded by 

n i E [ n r» n t] with *he no ™ na l value nj = .15 

no E [n^ ,nj] with the nominal value n® = .95 

€ [d 2 idy] with the nominal value d£ = .95 


P(3) := 

C(S) := 



Figure C. Interconnected Feedback System. 


Q( 3 ) 


Qx{ 3 ) _ 3 + 2 

Qii 3 ) 3+1 




AW 

AW 


P 11 (a) _ a 2 + 3 + 1 

Ai( s ) a 3 + 02 a 2 + 4a + ao 
P 12 (a) 6.6a 3 + 13.5a 2 + 15.5a + 20.4 

P 22 W a 3 + bjs 2 + 3.5a + 2.4 


and let the set of parameters p = [a 2 , ao, b 2 ] vary as 
follows: 


oj € ,41 w ith the nominal value o 2 = -3.0 

oo € [«o i°o] with the nominal value = 2.0 

&2 G [5 2 , ] with the nominal value & 2 = 3.5 

where [a 2 ,a 2 ] = [o 2 — w nj €, a% + u/ nj e] 

K.«o] =[aS- w »« e . a S + u; «oe] 

[6j ,&J] =[bj- w b , e, b° 2 + w bi e] 

with [w„ 2 = [1,1,1]. The Kharitonov ver- 
tex and segment polynomials as follows: 


E Pll (3,e) 

= 

K}Js t e) = K 3 P j3,e) 


= 

K* Pll (3,e) = P n (3,t) 

k 1 Pi2 M 

= 

K 2 p l3 (*,t) = K 3 Pl3 (3,t) 


— 

K* Pl3 (3,e) = P 12 (3 f e) 


— 



= 

% + 4a + aj a 2 + a 3 

K 3 P j3,e) 

= 

Kp 3l M 


— 

oj + 4a + oj a 2 + a 3 

Sp 2l W e ) 

= 

xhtM 

Sp 3l {3,e) 

— 

Sp 31 (a,e) 


= 

A^wo+a-A^c 

Sp 3l W e ) 

= 


KpJ s > e ) 

= 

k 2 P33 M 


= 

2.4 + 3.5a + 6ja 2 + a 3 

K 3 P33 (s,e) 

= 



= 

2.4 + 3.5a + b 2 a 2 + a 3 

Sp 33 ( a ’ € ) 

= 

Kp»M 

Sp 2 j( a ) e ) 

— 

Sp 33 {*,t) 


= 

A^„We) + (l-A)K 3 „(, 

Sp„W«) 

= 

K 3 pJ*,€). 


The sets of manifolds to be checked for stability are: 


Q 2 (s)[A a K> 31 (a, € ) + (1 - A 2 JJTjJa.e)] 
[A 2 K>„(a,e) + (l-A 2 )^„(a,e)]} 

We now solve the problem of checking these manifolds 
by overbounding this set by the convex hull of the 
vertices. After eliminating all duplicated segments, 
we have the following line segments which need to be 
checked. Note that the line segments listed below are 


functions of e and our objective is to find the max- 
imum value of e so that all E(s,c*) remain Hurwitz 
for allc G [0,c*]. 

E x {s t e) - {Qi{s)P 11 {s,€)P 12 {s i €) + 

<h(°)K 3 P Js,e)K 3 P j3,t), 

Q iWAi(*»0Aa(*.e) + 
Q2W*p„W «)*/•„('.«)) 
f?zW c ) = (Qi( a ) Ai(*» e ) Aj(*» c ) + 

Q2(*)K 3 p j3 t €)K%j3,€), 

Q 1 WA 1 WOA 2 WO + 
GaW*kW«)*ft,WO) 

E 3 {3,e) = (<?i(a)Pu(a,c)P 12 (s,e) + 

Q2(S)K 3 P J3,C)K 3 J3,C), 

Qi(a)Pn(a 1 e)Pi 2 (a 1 e) + 

Ei(s,() = (Qi(a)Pn(a,e)Pi 2 (a, e) + 

<?2W*P„W «)*£„(*»«)» 

Qi(a)Pn(a, e)Pi 2 (a, e) + 

A(a,e) = (Qi(a)Pu(a,e)P 12 (a,e) + 

Q 2 (a)^ Jl (a,c)K).„(a,€), 
Qi(a)Pn(a,e)Pi 2 (a,c) + 

Q 2 (a)K). Jl (a l£ )li:>„(a,e)) 
Ee(3,e) = (Q 1 (a)Pu(a,c)P 12 (a 1 €) + 

Q2{3)K 1 P j3, e )K 3 Pa ,(a,e), 
Qi(a)Pn(a,e)Pi 2 (a,e) + 

Q2(3)K 1 P j3,€)K 1 Pn (s,e)). 

Using a simple bisection method based algorithm we 
have 

e* = max{e|i^(5,e) 6 ft, Ve € [0,e*], i = 1, • • • , 6} 
= 0.6305. 

4. CONCLUDING REMARKS 

A simple technique to determine the parametric sta- 
bility margin for multilinearly dependent interval con- 
trol system has been given. The method is based 
on the multilinear generalization of the Box theorem. 
The method is powerful and reduces computational 
burden significantly. In fact, if used in conjunction 
with the Segment Lemma it completely eliminates 
frequency sweeping. Furthermore, by adding addi- 
tional vertices one may achieve an arbitrarily high 
level of accuracy. 
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Figure 1: Image Set (w =0.85) 



Figure 2: Convex Hull (a; = 0.85) 



Figure 3: Image Set (uj = 0.85) 



Figure 4: Convex Hull (u> = 0.85) 






